Analyzing fine scaling quantum effects on the buckling of axially-loaded carbon nanotubes based on the density functional theory and molecular mechanics method

In this paper, the quantum effects of fine scaling on the buckling behavior of carbon nanotubes (CNTs) under axial loading are investigated. Molecular mechanics and quantum mechanics are respectively utilized to study the buckling behavior and to obtain the molecular mechanics coefficients of fine-scale nanotubes. The results of buckling behavior of CNTs with different chiralities with finite and infinite dimensions are given, and a comparison study is presented on them. The differences between finite and infinite nanotubes reflect the quantum effects of fine scaling on the buckling behavior. In addition, the results show that the dimensional changes highly affect the mechanical properties and the buckling behavior of CNTs to certain dimensions. Moreover, dimensional changes have a significant effect on the critical buckling strain. Beside, in addition to the structure dimensions, the arrangement of structural and boundary atoms have a major influence on the buckling behavior.

Using MD simulations and continuum mechanics models, Yakobson et al. 9 studied the buckling of SWCNTs under axial compression.It was seen that the cylindrical shell configuration of a SWCNT with an aspect ratio of 6 would not change but axial compression makes it shorter.The effects of the axial compression and bending on the critical buckling strain of multi-walled carbon nanotubes (MWCNTs) has been studied by Chang et al. 58 using the molecular mechanics model.It was shown that only the part of the outer layers buckles first while the remaining inner part remains stable in a very thick MWCNT.In addition, investigation of the effects of tube size on the initial buckling wavelength revealed that the dependence of the initial buckling wavelength on the thickness of the MWCNT is low.Employing systematic MD simulations, Chang et al. 59 studied the bending of the thin (three-walled) and thick (seven-walled) MWCNTs at low temperature of 1 K.It was shown that the initial buckling mode of thick and thin MWCNTs are different from each other and only several outer layers buckle first, while the rest of the inner layers remain stable in a very thick MWCNT.Using comprehensive MD simulations, Jeong et al. 48,49 characterized the torsional behaviors of different CNT systems, such as bundled, chemically functionalized and nanotubes filled with other materials.Their results of Ref. 48demonstrated that the amount which affects the increase of the critical buckling moment relies on the number of inner tubes and the type of filling materials.The effect of initial stress on the torsional buckling behaviors of CNT systems was later examined by the authors 49 .They indicated that the critical torsional moment and stiffness could be significantly enhanced by the presence of initial stress.Using the nonlocal elasticity equations of Eringen, Pradhan and Reddy 60 investigated the buckling of SWCNT with Winkler foundation.It was shown that as the size of SWCNT decreases, the nonlocal effects increase and as nonlocal parameter increases, load ratio's decreases.In addition, further load ratio's increases with increase in Winkler modulus (K) for clamped-clamped, simply supported, clamped hinged boundary conditions.However, load ratio's decreases but with increase in Winkler modulus (K) for clamped free boundary condition.The buckling behavior of simply supported-simply supported single-walled carbon nanotubes (SWCNTs) was studied by employing the finite element method (FEM) by Bocko et al. 61 for the nanotubes with and without defects.The carbon nanotubes were modeled as beams and shells and it was seen that the critical buckling force of SWCNTs would decrease as the number of defects increases.It was also seen that the SWCNTs with the same diameter but different chirality have almost the same decrease in the critical buckling force.Using continuum mechanics models and molecular mechanics simulation, Ma et al. 62 studied buckling behaviours of the pre-stressed multi-walled carbon nanotubes (PS-MWCNTs) with two to six layers.They revealed three features of the buckling behaviour of PS-MWCNTs while the interlayer distance is considered as the key factor and the nanotubes are affected by the axial loading.It was demonstrated first, depending on the diameter of nanotubes, the buckling membrane force is not a monotonic function of interlayer distance.Second, the interlayer distance decreases for PS-MWCNTs with fixed intertube chirality as the buckling membrane force increases and third, the buckling membrane force increases as the number of walls increases for PS-MWCNTs with the same innermost tube.Moreover, they stated that the multi-shell continuum model and molecular mechanics simulation agree on the trend of the buckling membrane force as a function of interlayer distance, tube chirality index, and number of layers.Using molecular mechanics simulation surface Young's modulus of SWCNTs with different chiral angles and diameters were calculated by Fang et al. 13 .And based on same method Wan et al. 63 studied the mechanical properties of carbon nanotubes by developing a structural mechanics model.Their results stated that with the help of finite element analysis, the method they used is faster than atomistic simulations which makes it possible to simulate problems involving a large number of atoms.Using density functional theory calculations Chaudhuri et al. 64 investigated a series of zigzag and armchair nanotubes of carbon, boron and nitrogen with various values of tube diameters to understand the effect of the diameter values and the chirality on the energetics, structure and electronic properties of nanotubes.It was revealed that carbon-boron (CBNT) and carbon-nitrogen (CNNT) nanotubes mainly shows metallic behavior based on the composition.While boron-nitrogen (BNNT) nanotubes seemed to have semiconducting behavior.In addition, the stability of the nanotubes were proven to be dependent on the respective chiralities.Hwang et al. 65 investigated the mechanical behavior of multi-walled carbon nanotubes (MWNTs) utilizing the molecular dynamics (MD), under uniaxial tensile loading while MWNTs are fixed at both ends.The range of Young's modulus of these nanotubes was found to be in the range between 0.85 and 1. 16 TPa.Different radii and lengths of the tubes were stated as a possible reason to the difference in the obtained modulus.It was also found that changing the number of boundary layers would have no effect on the calculated mechanical properties, which indicates indicating the role of mechanical boundary conditions in the MD simulations.Cao 66 developed carbon fiber-reinforced thermoplastic sandwich composites and introduced an innovative technique to enhance the interphase between the skin and core.This method involved the incorporation of resin-coated CNT-yarn fillers, offering a viable solution to improve the debonding toughness in low-weight-gain sandwich composites.The examination of fracture characteristics in reinforced concrete incorporating carbon nanotubes (CNTs) was conducted through the utilization of double-K fracture parameters by Zheng et al. 67 .The findings indicated that the incorporation of CNTs can enhance the material's resistance to cracking.Kim et al. 68 systematically investigated alterations in the fiber microstructures, with a particular focus on the inter-bundle and intra-bundle voids during the process of solution spinning.Near the spinneret exit, it becomes evident that extensional deformation, achieved through drawing, plays a crucial role in orienting the fibers.Monavari et al. 69 investigated the electronic response of single-walled carbon nanotubes (SWCNTs) and a carbon nanobelt (CNB) to N-linked and O-linked SARS-CoV-2 spike glycoproteins, using ab initio quantum mechanical approach.The CNTs were selected from three zigzag, armchair, and chiral groups.They examine the effect of carbon nanotube (CNT) chirality on the interaction between CNTs and glycoproteins.Results indicated that the chiral semiconductor CNTs clearly response to the presence of the glycoproteins by changing the electronic band gaps and electron density of states (DOS).Since the changes in the CNTs band gaps in the presence of N-linked were about two

DFT simulation details
The current study is based on the work of Ansari et al. 8 , which explained how the final form of equations used in this paper are computed.Thus, here we just present the necessary equations to avoid duplications.The mechanical properties of the graphene sheet based on the strain energy calculation is Acquired through Quantum-Espresso coding procedure 73 .The calculations are based on the density functional theory [74][75][76][77][78][79][80][81][82][83][84][85][86][87][88][89] and generalized-gradient approximation (GGA) function is applied along with the Perdew-Burke-Ernzerhof (PBE) exchange correlation 90,91 and the projector-augmented wave (PAW) potential is applied were employed for the self-consistent total energy calculations and geometry optimizations.After optimizing the initial parameters of structures' input, the kinetic energy cutoff for the plane wave basis set has been converged to 80Ry .In order to compute kinetic energy cutoff, total energy variations versus cutoff energy variations have been computed and we plot the variations.Finally, in cutoff energy 80Ry , difference in total energy versus cutoff energy reached lower than 0.005 eV.Hence, we selected 80Ry as our optimized cutoff energy.The same trend we employed on optimizing k-point grid.The Brillion Zone was sampled using a 20 × 20 × 1 Monkhorst Pack k-point grid.The basic concepts and relations used in the Quantum-Espresso code are expressed in the literature 8 .

Potential energy
The total potential energy,E t , in the empirical force field method of molecular mechanics, can be stated as the sum of several energies due to valence of bonded interactions or bonded and non-bonded interactions (Fig. 1): In which U ρ , U θ , U ω , and U τ are energies associated with bond stretching, bond angle variation, bond inversion, and torsion, respectively; U vdw and U es are also associated with van der Waals and electrostatic interactions, respectively 11,[92][93][94][95] .Depending on the particular materials and loading conditions considered, different functional forms might be utilized for these energy terms.For the present buckling problem of carbon nanotubes, the (1) terms U ρ , U θ and U ω are expected to be significant in the total potential energy of the system.Furthermore U vdw can be neglected because SWCNTs are considered herein.As Hooke's law has proven to be efficient and precise enough to describe the behavior of atoms under small deformation 94 , it is frequently employed to characterize the interactions between bound atoms in the system, therefore, using Ref. 96 , Eq. (1) may be written in the form: Here the subscript k = 1, 2, 3 .Besides, it should be remarked that the constant coefficient 1 2 of the first term in Eq. ( 2) is to ensure that bond stretching energy is considered only once.

Axisymmetric buckling of single-walled chiral nanotubes
The structure of a single-walled carbon nanotube is often characterized by a pair of integers ) n, m (, representing its helicity 97 .The geometrical parameters are same as Sect."Molecular mechanics model" of Ref. 96 .A carbon nanotube under compression tends to buckle when its axial strain exceeds the critical compressive strain,ε 0 .Prior to buckling, the variation of bond length, bond angle and inversion angle from equilibrium values are obtained as Eqs.( 9) to (27) of Ref. 8 .Then, using Eq. ( 24) of the mentioned reference, one would have: dr ij2 is also obtained as same as above approach: dr ij3 is obtained By means of Eqs. ( 13)-( 15), ( 16), ( 18), (21) of Ref. 8 : (2) (3) www.nature.com/scientificreports/According to the displacement of the atoms of SWCNT after buckling, the variations of bond length ( r ij1 , r ij2 , r ij3 ), bond angle ( �θ ij1 , �θ ij2 , �θ ij3 ), and inversion angle ( �β ij ) of the deformed SWCNT after buckling can be obtained.These quantities are given in the Appendix A of Ref. 96 . in which the following parameters are used, where R , the tube radius can be expressed in the form, where ξ i,j the radial displacement of the atom ij , is assumed to be small in comparison with r 0 .For atoms located at type B position as shown in Fig. 5 of the Ref. 96 , the term ξ i,j+1 in the preceding equations must be replaced by ξ i,j−1 .In fact, there is no crucial difference between atoms at type A and type B positions.For convenience, the mathematical procedure described herein is only for atoms at a type A position.

Stability equation
Making use of the parameters discussed in Sect."DFT simulation details", the potential energy of the system prior to buckling, dE , and after buckling, E , can be computed with using Eq. ( 2) The free energy of the system, , is, The condition for to be extremum requires,

Buckling strain
In the case of axisymmetric buckling, the radial displacement of atom ij is of the form In which L is the tube length, x ij is the longitudinal coordinate of atom ij , and Z 0 , the rigid radial extension of the nanotube under a compressive strain of ε 0 12 .Introducing Eq. ( 12) into (13) yields the stability equation from which the buckling strain can be attained for its nontrivial solutions, Z = 0.
Force constants ( K ρ , C θ and C ω ), Young's modulus, Poisson's ratio and flexural rigidity of fine-scale structures To investigate the buckling behavior and other properties related to the structural energy of materials using molecular mechanics, the most important factor is to derive the coefficients of these equations in order to find the changes in energy values due to changes in material size and structural atomic arrangement.In this section, in order to account the fine-scale nature of the materials on its structural properties and to study buckling behavior and other energy related issues that can be extracted from molecular mechanics equations, the coefficients of appropriate molecular mechanics equations must first be extracted.To achieve this, due to the properties of fine-scale nanosheets, which leads to the creation of nanotubes with specified length, the molecular mechanics (7) .
coefficients related to fine-scale structures have been extracted.The characterization of graphene as a twodimensional carbon nanostructure that can be used as a standard for comparing the change of properties with respect to dimensional changes has also been studied.Structural energy variations were nonlinear and more varied widths should be selected in areas where changes are more frequent and in areas where energies are closer to 2D structures, points should be selected in accordance with this trend.For this purpose, we first need to select specific widths of the nanosheets along the desired directions, which can be used to characterize the properties change process as well as to optimize the calculations.As a result, in order to achieve this purpose, there are some widths considered in this paper whose shape are presented in Fig. 2 and their properties are presented Table 1.Widths 2a, 3a, 5a and 8a are selected in the armchair direction and 2z, 3z, 5z and 9z in the zigzag direction so that the widths are as close as possible to each other in terms of structure characteristic so that a closer comparison can be made between the effects of different directions.It should be noted that the points were chosen so that they would have higher abundance in less than 20 angstroms, where the energy variations were with respect to the size difference noticeable.Also around and higher than 20 angstroms, some points are selected to compare the properties and extent of changes.After determining the mechanical properties of the plate nanostructures according to the references 8,96 , the molecular mechanics coefficients are obtained, and the results of the mechanical properties and molecular mechanics coefficients are presented in Table 2.

Numerical results and discussion
Using the following formula buckling mode parameter ξ can be obtained: www.nature.com/scientificreports/Here, L, r 0 and m 1 represent the tube length, C − C bond length and the number of half waves in the axial direction respectively.It is clear that the minimum value of the buckling strain (i.e. the critical buckling strain) is sensitive to the tube diameter so that the critical buckling strain decreases with the increase of tube diameter.By introducing the buckling wavelength It can be observed that the critical buckling strain occurs at higher buckling mode parameter (or at lower wavelength), as previously reported in Ref. 12 .
In Fig. 3, the values of the buckling strain (ε 0 ) of infinite tubes with zigzag atomic arrangement and chirality of (n, 0) are shown in terms of inverse mode modulus (ξ ) .The results show that by varying the buckling mode parameter from maximum value to minimum, the buckling strain values also change to a minimum value.This value can be considered as the critical buckling strain.The results also show that for the larger buckling mode parameters, the maximum buckling strain values are larger than the similar values for the smaller buckling mode parameters over the study period.In addition, for nanotubes with zigzag atomic arrangement, the buckling strain values change from a smaller and closer range to a wider range with decreasing buckling mode parameter value.It can be seen that as the diameter of the nanotube structure increases, the extent of buckling strain increases and the buckling strain values decrease.In smaller buckling mode parameters, the buckling strain is more dependent on the diameter and decreases with increasing the diameter.Moreover, it is seen that by increasing the diameter of nanotube, critical buckling strain occurs in smaller buckling mode parameter and the critical buckling strain decreases.In addition, as the diameter of the buckling strain changes, the buckling mode becomes closer to each other.Furthermore, it is seen that in the larger buckling mode parameters the variation of changes in the buckling strain are higher and lower in the smaller buckling mode parameters.
The results of variations in buckling strain values of zigzag nanotubes with specified lengths obtained by wrapping armchair nanosheets are presented in terms of inverse buckling mode parameter are plotted in Fig. 4. By comparing the results, the buckling strain susceptibility of zigzag nanotubes from dimensional changes and quantum effects of fine scaling will be determined.The results show that by increasing the length of the structure in a well-defined process, the buckling strain values start from smaller values and increase as the length of the nanotubes increase.For smaller nanotubes, the buckling strain values for the larger buckling mode parameters have a wider range of values, which increase by increasing buckling mode parameter.It can be concluded that as the length of the structure increases, the range of the buckling strain values increases.In other words, (15)   by increasing the length of the structure and the diameter of the structure, the buckling strain is affected by the quantum effects of fine scaling, which is greater for the diameter variations compared to the longitudinal variations.The curvature effects can also be considered effective in this case.The results show that due to variations in the buckling mode parameter, the range of buckling strain values in nanotubes with larger lengths have increasing trend.The results show that for the nanotubes with specified length, the values of buckling strain and buckling mode parameter, experiences a decreasing trend with higher rate of changes with increasing the diameter of nanotube, which is same as the nanotubes with infinite length.As the length of structure increases the rate of change also increases and changes in the values of buckling strain get closer to each other at larger www.nature.com/scientificreports/diameters.Furthermore, it can be seen that as the diameter of the structure increases, the changes in the values of buckling strain experiences a decreasing trend as the buckling mode parameter changes.In Fig. 5, the variations in the buckling strain values of zigzag nanotubes with specified lengths obtained by wrapping zigzag nanosheets are shown in terms of the inverse buckling mode parameter.The results show that, same as the nanotubes with specified lengths obtained by wrapping armchair nanosheets and nanotubes with infinite length, by changing the buckling mode parameter, the buckling strain has a minimum value.In addition, it is seen that as the diameter of the structure increases, the variations in the buckling strain values become larger, and the changes in the buckling strain values get closer to each other for different nanotubes as the diameter of the nanotubes increases.In other words, the trend of changes in the buckling strain of nanotubes is similar to the buckling mode parameter changes in nanotubes with specified lengths and the ones with infinite lengths, but the quantum effects of fine scaling would result in of these values.Moreover, the results also show that by increasing the length of structure for the nanotubes obtained by wrapping zigzag nanosheets, the buckling strain increases and the magnitude of changes is higher in smaller diameters and compared to the nanotubes obtained by wrapped armchair nanosheets, the magnitude of the changes is larger.Furthermore, by increasing the length and diameter of the structure, the trend of change in the buckling strain with respect to the buckling mode parameter also decreases.
The buckling strain values (ε 0 ) with respect to the inverse mode modulus (ξ ) are plotted in Fig. 6 for nanotubes with infinite length, armchair atomic arrangement, and chirality of (n, n) .The results show that in the domain of the buckling mode parameter from small to large, the buckling strain values for the nanotubes with infinite length and armchair atomic arrangement, changes toward having a minimum value, which is the critical buckling strain.The smaller buckling mode parameter has a larger range of change in the buckling strain values and as the buckling mode parameter increases, changes occur in a smaller range.In addition, as the diameter of the nanotubes increases, the buckling strain values decrease and these values become closer for the nanotubes with larger diameters.The trend of buckling strain changes for nanotubes with infinite length and armchair atomic arrangement is similar to that of zigzag nanotubes, except that the buckling strain values of armchair nanotubes are in smaller domain with higher rate of change compared to zigzag nanotubes.
Figure 7 shows the variations in the buckling strain values with respect of the inverse mode buckling for armchair nanotubes with the specified lengths obtained by wrapping armchair nanosheets.The results show that by changing the buckling mode parameter, the nanotube buckling strain starts with changes in values close www.nature.com/scientificreports/ to each other in the larger buckling mode parameters and then when the buckling mode parameter decreases, it shows a wider range of changes which is similar to armchair nanotubes with infinite length.However, the quantum effects of the fine scaling and shape of the nanotubes have results in reduced buckling strain values as the length of the nanotube structure decreases.The range of buckling strain variations for different diameters in the larger buckling mode parameters varies over a larger range and as the length of the structure increases, the buckling strain values become closer to each other.It can be seen that as the length and diameter of the nanotubes increase, the buckling strain for a specific diameter, changes with lower slope.By comparing the changes in the quantities of the fine-scale armchair and zigzag nanotubes, it should be noted that the overall trend is similar.
Except that the critical strain values of the armchair nanotubes are lower than the zigzag nanotubes and the www.nature.com/scientificreports/quantum effects caused by fine scaling, is higher in nanotubes with armchair atomic arrangement compared to the ones with zigzag atomic arrangement.The variations in the buckling strain values with respect to inversed buckling mode parameter are plotted in Fig. 8 for the armchair nanotubes with specified lengths obtained by wrapping zigzag nanosheets.The results show that as the length of the structure increases, the buckling strain has an increasing trend and these changes are greater for smaller lengths.According to the results, similar to infinite armchair nanotubes, as the buckling mode parameter changes from larger values to lower values, the buckling strain changes in a similar trend as it starts with values with less variation range and eventually the range of variation becomes wider.In addition, as the diameter of the nanotubes becomes larger in higher buckling mode parameters, the buckling strain values become closer to each other and as the parameter of buckling mode becomes smaller, the buckling strain values would have a wider range of variation.Compared to infinite armchair nanotubes, the buckling strain of armchair nanotubes obtained by wrapping zigzag nanosheets, has smaller values in the larger buckling mode parameters, and with the variations in the diameter and the buckling mode parameter, the buckling strain values have more variability.The graph also shows that by decreasing the buckling mode parameter, the variations of the buckling strain values are such that the larger buckling strain values are more on the side of the parameters with the smaller buckling mode.Compared to armchair nanotubes obtained by wrapping armchair nanosheets, the ones that are obtained by wrapping zigzag nanosheets have smaller buckling strain values, greater buckling strain values with respect to changes in diameter and larger buckling mode parameter.Compared to zigzag nanotubes obtained by wrapping zigzag nanosheets, it also has smaller buckling strain values, greater buckling strain values with respect to diameter and larger buckling mode parameter variations, and a larger slope of buckling strain values in smaller buckling mode parameters.The results show that with increasing diameter and length of structure, the changes of buckling strain with respect to buckling mode parameter have a decreasing trend.
In Fig. 9, the variations in the buckling strain values of the infinite nanotubes with different chirality are plotted with respect to the inverse buckling mode parameter.It can be seen that the trend of changes in the buckling strain values is in a way that with changes in the larger buckling mode parameter, the buckling strain values start from values within a close range, and a similar trend and eventually in the larger buckling mode parameters they would have a wider range of values.In addition, as the diameter increases, there is also an increase in the range of the buckling strain values in the smaller buckling mode parameters.Moreover, it is found that the buckling strain values of infinite chiral nanotubes in the larger buckling mode parameters have smaller buckling strain values compared to the infinite zigzag nanotubes.In addition, the range of buckling www.nature.com/scientificreports/strain values is such that in some smaller diameters, the buckling strain values are larger and in the larger diameters, the buckling strain values are lower, compared to armchair nanotubes.In the smaller buckling mode parameters, the chiral nanotubes with smaller diameters have wider range of values and in larger diameters it is reversed compared to the armchair and zigzag nanotubes.The results show that as the diameter of the structure increases, the slope of the buckling strain changes with respect to the buckling mode parameter experiences a decreasing trend.The trend of changes in the buckling strain values of chiral nanotubes, compared to armchair and zigzag nanotubes, are larger, in the higher buckling mode parameters, and smaller in the lower buckling mode parameters.Due to the different chirality of the these nanotubes, the changes of their buckling strain values with respect to the diameter changes are less regular compared to the zigzag and armchair nanotubes, but the trend of changes in buckling strain values with respect to the buckling mode parameter is similar.
The variations in the buckling strain values with respect to inverse buckling mode parameter are plotted in Fig. 10 for nanotubes with different chirality and specified lengths obtained by wrapping armchair nanosheets.The results show that with increasing the length of the structure, the buckling strain values experiences an increasing trend, and this trend is reduced with increasing the length.Compared to chiral nanotubes with infinite length, the ones with specified length have smaller buckling strain values in larger buckling mode parameters that fall within the range of infinite nanotubes as the length of nanotubes increases.Moreover, the results show that as the diameter and length of the structure increase, the variation of the buckling strain with respect to the buckling mode parameter experiences a decreasing trend.In addition, the slope of graph in finite length nanotubes is higher in lower buckling mode parameters compared to nanotubes with infinite length.Minimum buckling strain values of finite length nanotubes increase from larger buckling mode parameters to smaller parameters as the length increases.Compared to Armchair and Zigzag nanotubes with finite length, the results show that the buckling strain values of different chiral nanotubes are smaller than the buckling strain values in larger buckling mode parameters for the finite length zigzag nanotubes and it is close to the values of armchair nanotubes with finite length.This trend continues at smaller buckling mode parameters except that at smaller buckling mode parameters the trend of change in the buckling strain values of the larger chiral nanotubes is higher.
Figure 11 shows the variations in the buckling strain values with respect to the inverse of the buckling mode parameter for nanotubes with different chirality and specified lengths obtained by wrapping zigzag nanosheets.The results show that with increasing the length of the structure, the buckling strain values have an increasing trend and these changes decrease with increasing the length.Compared to infinite chiral nanotubes, chiral nanotubes with specified length have smaller buckling strain values in the larger buckling mode parameters, and compared to chiral nanotubes with specified length obtained by wrapping armchair nanosheets, the buckling strain values in the larger buckling mode parameters are in close range to each other.Compared to zigzag and armchair nanotubes with specified lengths obtained by wrapping zigzag nanosheets, the buckling strain values at the larger buckling mode parameters were close to the buckling strain values of the armchair nanotubes and smaller than the zigzag nanotubes obtained by wrapping zigzag nanosheets.Due to the chirality and the way the atoms are positioned in the structure of the chiral nanotubes, by varying the buckling mode parameter and the diameter of the nanotubes, their buckling strain values their buckling strain values are less regular compared to other nanotubes.Moreover, compared to other nanotubes, the buckling strain values in the smaller buckling mode parameters have a higher slope and the minimum buckling strain values are placed on the side of the smaller buckling mode parameters.Furthermore, as the diameter and length of the structure increase, the changes in the buckling strain with respect to the buckling mode parameter have a decreasing trend.
In Fig. 12, the results of critical buckling strain analysis with respect to diameter of nanotube are presented for nanotubes with infinite length and different chiralites.Here, diameters up to 30 angstroms have been studied to investigate the quantum effects of fine scaling on the critical buckling strain of nanotubes.The results show that zigzag nanotubes with infinite length in the range of 5 to 30 angstroms have the highest critical buckling strain compared to nanotube with other chiralites.Critical buckling strain variations of zigzag nanotubes with infinite length have a uniform trend of change within the specified range for diameters due to their structural www.nature.com/scientificreports/features and atomic arrangement.This uniform trend of change means that as the diameter of the nanotubes increases, their buckling strain values decrease and this trend is in a way that at diameters with the values less than 15 angstroms, larger changes would occur in the buckling strain compared to those with diameter values greater than 15 angstrom.In other words, as the diameter of the zigzag nanotubes increases, the buckling strain continues to decrease but the amount of these changes decreases, indicating that increasing the diameter reduces the resistance of the structure to the axial loading.In addition, at larger diameters where the quantum effects of fine scaling decrease, the trend of changes in the buckling strain also decreases, so some of the changes that lead to the reduction of the buckling strain can be attributed to the quantum effects of fine scaling.Moreover, the results show that the trend of change in the buckling strain of zigzag nanotubes with infinite length is always decreasing at the small diameters.For the armchair nanotubes with infinite lengths within the range of 5 to 30 angstroms, it can be seen that the critical buckling strain also shows a decreasing trend, except that the buckling strain of infinite armchair nanotubes in the range of 5 to 10 angstroms has different trend of changes compared to other diameters.In this range, the results of the critical buckling strain is closer to the results of zigzag nanotubes.In the range of 10 to 15 angstroms, the trend of changes for the critical buckling strain is different from the other ranges and the results are more distant from the critical buckling strain of the zigzag nanotubes.However, the results are closer to the critical buckling strain of nanotubes with chirality of (n, n/2) than the zigzag nanotubes.The trend of change in the critical buckling strain results of the armchair nanotubes will change once more for diameters higher than 15 angstroms.As the diameter increases, the critical buckling strain results moves away from the results of zigzag nanotubes and gets more closer to nanotubes with chirality of (n, n/2) and eventually from diameters higher than 25 angstroms, it is quite close to the results of nanotubes with chirality of (n, n/2) .
The results also investigate the nanotubes with chirality of (n, n/2) in the diameter range of 5 to 30 angstroms.It can be seen that the critical buckling strain of nanotubes with chirality of (n, n/2) has a decreasing trend in the studied area and with increasing diameter, the reduction intensity of the critical buckling strain decreases.The results show that for diameters in the range of 5 to 10 angstroms, the buckling strain values of nanotubes with chirality of (n, n/2) are closer to the results of armchair and zigzag nanotubes, and at larger diameters, it differs from the results of the armchair and zigzag nanotubes and the maximum spacing from zigzag nanotubes is reached in 30 angstroms.However, at the diameters larger than 25 angstroms, the results are close to the results of the armchair nanotubes.In general, as the diagram shows, the critical buckling strain of nanotubes with different chirality are obtained between the results of zigzag nanotubes and nanotubes with chirality of (n, n/2) .
The results show that at very small diameters the buckling strain values of the nanotubes with different chirality  are closer to each other.This may be due to the very small diameters of the nanotubes, and the type of buckling that takes place.Then, at slightly larger diameters, the trend is such that the buckling strain of nanotubes with different chirality differs markedly, and this area is where the chirality and arrangement of structural atoms have a greater effect on the buckling of nanotubes with different chirality.Subsequently, at larger diameters, that the structure of nanotubes slowly become closer to the structure of their constituent sheet, with slight increase in the diameter the critical buckling strain of the nanotubes with different chiralities become close to each other and the buckling occurs closer to the buckling of the sheet.www.nature.com/scientificreports/ Figure 13 shows the variation of critical buckling strain versus carbon nanotube diameter for the state of armchair.Also obtained results from Chang's molecular mechanic 12 and Yakobson's continuum method 9,20 are presented in Fig. 13 which shows large difference with those of molecular dynamic results.It can be seen that Chang's model 12 tends to overestimate the buckling strains, especially when the tube diameter decreases.In contrast, the results obtained from the present analysis are found to be in excellent agreement with the ones from the molecular dynamic simulations, which indicates the capability of the present approach in predicting buckling strains of carbon nanotubes.
Table 3 is provided to show the critical buckling strains of single-walled CNTs as obtained by various researchers using atomistic simulations under different simulation conditions.It's clear that the obtained results in this study is in good agreement with those obtained by molecular dynamic method in comparison with continuum method and molecular mechanics done ever before 9,10,51,57,[98][99][100][101][102][103][104][105][106] .Yakobson et al. 9 obtained a critical buckling strain of 5% for a (7, 7) tube from classical molecular dynamics (MD).This value is much lower than Chang's prediction of about 10.7% from molecular mechanics 12 and the prediction by Yakobson's continuum model of 7.7% 9,20 and what obtained in present prediction for critical buckling strain is 6.1% that has less difference with Yakobson's results in comparison with Yakobson's continuum model and Chang's molecular mechanics.For (10,10) tube the Chang's prediction is 6.92% and Yakobson's continuum model prediction is 5.69% which is more than Zhang's 104 results obtained by molecular dynamics (MD) which is 4.29% and the obtained critical buckling strain in this prediction is 4.23%.Differences between Chang's molecular mechanics and present molecular mechanics method are due to technique of acquisition of energy equation's terms.The present approach gives the total potential energy without any simplification.
To investigate the influence of the chirality on the buckling strains of nanotubes, the variation of critical buckling strain versus diameter of nanotube is plotted in Fig. 13 for armchair, zigzag and chiral tubes.The difference between results obtained for armchair and zigzag tubes becomes more pronounced when the diameter of tube increases.As it can be seen in this figure, the stability of zigzag carbon nanotube under axial load is better than armchair carbon nanotube and the minimum stability is for (n, n/2) tube.Furthermore, this figure presents the critical buckling strain of chiral nanotube in which the variation of critical buckling strain in chiral tube is such that for a certain diameter, increasing in the chiral angle from θ = 0 to θ = π/12 , critical buckling strain changes from zigzag to (n, n/2) tube.With continuing the trend of increase in chiral angle from θ = π/12 to θ = π/6 , critical buckling strain changes from ( n, n/2 ) tube to ( n, n ) tube.
Figure 14 shows the results of the critical buckling strain with respect to nanotube diameter for finite length nanotubes with different chiralities obtained by wrapping armchair nanosheets.The critical buckling strain variations of finite length nanotubes are in a trend similar to the critical buckling strain of nanotubes with infinite length.As the results show, in general, the changes are in a way that the critical buckling strain decreases with increasing the diameter of the structure, which depends on the length of structure, diameter and the arrangement of the atoms.In all nanotubes of different lengths, the critical buckling strain of the zigzag nanotubes has a decreasing trend with respect to the diameter, and with increasing structure length, the critical buckling strain of the zigzag nanotubes has a decreasing trend compared to the shorter lengths.By increasing the length of zigzag nanotubes, the trend of change in their critical buckling strain decreases, so that for shorter lengths and for diameters of lower than 15 angstroms, the critical buckling strain variations are greater.For armchair nanotubes, www.nature.com/scientificreports/as the length of the structure decreases, the critical buckling strain decreases and by increasing the length of the structure, the values of the critical buckling strain decrease at a specified diameter.At shorter lengths, the critical buckling strain of armchair nanotubes in the larger diameter range is closer to the critical buckling strain of the zigzag nanotubes and by increasing the length of the nanotubes at smaller diameters, it deviates from the critical buckling strain of the zigzag nanotubes.The critical buckling strain of the armchair nanotubes at smaller lengths, even if it has a large diameter, is still smaller than the critical buckling strain of nanotubes with the chirality of (n, n/2) .However, as nanotubes grow larger in diameter, they deviate from the results of zigzag nanotubes and become closer to the results of (n, n/2) nanotubes.The buckling strain of the armchair nanotubes in smaller lengths has uniform decreasing trend and as the length increases, the intensity of this trend and the critical buckling strain changes would be different from the variations in smaller lengths.It can be seen that for (n, n/2) nanotubes the critical buckling strain has a similar trend to the infinite nanotubes and compared to the ones with infinite length, the trend of changes and the decrease in critical strain values reduces as the length increases.Moreover, the critical strain values for a specified diameter experiences a decreasing trend as the length of nanotube increases.Comparison of the critical buckling strain results of the finite and infinite length nanotubes reveals that the as the length the diameter of the nanotubes become smaller, the critical strain values of the nanotubes with different chirality become closer to each other.In addition, at larger lengths and larger diameters, chiral differentiation becomes more evident on the trend of changes of critical buckling strain, insofar as increasing the diameter does not eliminate the curvature effects of the nanotubes.In Fig. 15, results of the critical buckling strain with respect to nanotube diameter of finite-length nanotubes obtained by wrapping zigzag nanosheets are plotted.The results show that the general trend of buckling strain variations of nanotubes with different chiralities obtained by wrapping zigzag nanosheets is similar to finite and infinite nanotubes obtained by wrapping armchair nanosheets in a way that with increasing the diameter of the www.nature.com/scientificreports/structure, the critical buckling strain decreases.However, the variations in their critical buckling strain values are different.The buckling strain variations of finite zigzag nanotubes obtained by wrapping zigzag nanosheets initially begin with values lower than the critical buckling strain values of zigzag nanotubes with infinite length and zigzag nanotubes obtained by wrapping armchair nanosheets and as the diameter increases, it experiences a decreasing trend.In addition, with increasing length, the buckling strain values increases in specific diameter until it reach the values of infinite nanotubes.The slope of the critical buckling strain variations at smaller diameters is greater than that of the larger diameters, and as the diameter changes, there is a regular decreasing trend in the values.The trend of changes in the critical buckling strain values of zigzag nanotubes obtained by wrapping zigzag nanosheets is lower than that of infinite and finite zigzag nanotubes obtained by wrapping armchair nanosheets.The results show that the critical buckling strain of nanotubes obtained by wrapping zigzag nanosheets changes with increasing diameter in a way that, at shorter lengths, the values are smaller than that of infinite nanotubes and armchair nanotubes obtained by wrapping armchair nanosheets with similar diameters.In addition, with increasing length, it approaches critical buckling strain values of nanotubes with infinite length.
As the diameter increases, the buckling strain of the armchair nanotubes are initially closer to the critical buckling strain values of the zigzag nanotubes and as the diameter increases, the changes are such that they approach the critical strain of(n, n/2) nanotubes.Moreover, the results show that the critical buckling strain of (n, n/2) nanotubes obtained by wrapping zigzag nanosheets, starts with a similar trend of change and lower values than infinite (n, n/2) nanotubes obtained by wrapping armchair nanosheets which decreases with increasing diameter and with increasing length, it continius to enlarge until it reaches the value of infinite nanotubes.
In Table 4, the critical buckling strain results of zigzag nanotubes with finite and infinite lengths obtained by wrapping armchair and zigzag nanosheets are presented.In order to study and quantify the critical buckling strain changes of the nanotubes, the percentage of their critical buckling strain changes has also been studied.To compare the results, zigzag nanotubes with diameter from about 2 to 30 angstroms were studied.The results show that in zigzag nanotubes with infinite length, critical buckling strain decreases with increasing diameter, such that in smaller diameters the trend of changes is higher than the changes in larger diameters and in larger diameters the trend of change in the critical buckling strain is reduced due to the reduction of nanotube curvature.The results show that at the diameters less than 6 angstroms, the percentage of buckling strain changes is greater than the percentage of diameter changes, which then decreases.The results show that the ratio of critical buckling strain variations to diameter changes in diameters larger than 25 angstroms is less than one.This implies that in finite length nanotubes, the buckling strain is more dependent on the diameter changes and because of the small length of the structure, the quantum effects caused by the changes in the diameter to the larger diameters are even more significant.Results for zigzag nanotubes obtained by wrapping 3a armchair nanosheets show that the buckling strain variations compared to nanotubes with shorter length are lower, and these variations are higher in larger nanotubes with larger lengths.The results also show that the ratio of critical buckling strain to nanotube diameter changes at diameters below 15 angstroms is higher than one, which decreases for larger www.nature.com/scientificreports/diameters.Moreover, the results show that with increasing length of zigzag nanotubes obtained by wrapping armchair nanosheets, the critical buckling strain variations at smaller diameters decrease and it gets closer to the variations of nanotubes with infinite lengths.In addition, as the diameter increases, trend of changes of the critical buckling strain reduces and eventually approaches the trend of changes of the critical buckling strain of the nanotubes with infinite length.The ratio of the critical buckling strain changes to nanotube diameter changes in zigzag nanotubes obtained by wrapping 5a and 8a armchair nanosheets at diameters less than 10 angstroms is higher than one and it is reduced at larger diameters.For zigzag nanotubes obtained by wrapping 2z zigzag nanosheets, the results show that the critical buckling strain at larger diameters has more variations compared to diameter changes and with increasing the diameter the trend of changes decreases.Compared to infinite length nanotubes, the critical buckling strain variations are greater and compared to the finite length nanotubes obtained by wrapping 2a armchair nanosheets at smaller diameters they have smaller buckling strain variations, which is reversed at larger diameters.The ratio of critical buckling strain variations to nanotube diameters, in diameter less than 25 angstroms is higher than one which decreases in larger diameters.For zigzag nanotubes obtained by wrapping 3z zigzag nanosheets, the buckling strain variations are smaller than those of with shorter length and nanotubes obtained by wrapping 3a armchair nanosheets, and are closer to the variation of infinite length nanotubes.Furthermore, the results show that the critical buckling strain variations of zigzag nanotubes obtained by wrapping zigzag nanosheets are lower than zigzag nanotubes obtained by wrapping armchair nanosheets at smaller diameters, which increases at larger diameters.In addition, with increasing the length of the nanotube structure, due to the nature of the sheets which are wrapped to form theses nanotubes, the critical buckling strain changes occur from lower to higher in smaller diameters.By increasing the length of zigzag nanotubes obtained by wrapping zigzag nanosheets, the trend of changes of the critical buckling strain reduces with respect to the changes of the diameter.In addition, the ratio of the critical buckling strain to the diameter changes of the nanotubes in the zigzag nanotubes obtained by wrapping 5z and 9z zigzag nanosheets is larger than one in diameters lower than 10 angstroms, which reduces in larger diameters.The critical buckling strain results of armchair nanotubes with infinite and finite lengths obtained by wrapping armchair and zigzag nanotubes are presented in Table 5.The results show that with increasing the diameter of nanotubes, the critical buckling strain has a decreasing trend.For this study the diameters in the range of about 2 to 30 angstroms are chosen for armchair nanotubes.In the infinite armchair nanotubes, with increasing the diameter, the critical buckling strain reduces and as the diameter increases, the trend of changes also increases in a way that the ratio of critical buckling strain to nanotube diameter changes has always been greater than one and increases until the diameter reaches the value of 25 angstroms.Compared to zigzag nanotubes with infinite length, the critical buckling strain variations were higher for the armchair nanotubes, and this increase in the magnitude of the variations with respect to changes in the diameter continues up to larger diameters.For the finite length armchair nanotubes obtained by wrapping 2a armchair nanosheets, the results show that with increasing diameter the critical buckling strain structure of the nanotubes has a decreasing trend.In addition, compared to the infinite length armchair nanotubes, it can be seen that in the critical buckling strain starts at larger values in the smaller nanotubes and in the larger diameters the critical buckling strain of the smaller length nanotubes is greater.Then, with increasing the diameter, for diameters larger than 15 angstroms, the critical buckling strain variations of the nanotubes with infinite length are increased.Compared to similar zigzag nanotubes with finite length, the critical buckling strain with respect to similar diameters have lower values and would experience greater variations with respect to diameter changes.The ratio of critical buckling strain to nanotube diameter changes has always been greater than one, but with increasing the diameter, unlike zigzag, armchair, and infinite length nanotubes, it experiences a decreasing trend, which is reversed at diameters greater than 20 angstroms.For limited length armchair nanotubes obtained by wrapping 3a armchair nanosheets the critical buckling strain has a decreasing trend with increasing diameter, and the buckling strain values are higher than those of nanotubes with infinite length and lower than the smaller length nanotubes.Critical buckling strain changes have an increasing trend and it has fewer changes compared to nanotubes with shorter length.The ratio of critical buckling strain to nanotube diameter changes has always been greater than one and in the diameters less than 10 angstroms, it has a decreasing trend.Compared to zigzag nanotubes made from similar nanosheets, they have smaller buckling strain and more variations.For armchair nanotubes obtained by wrapping 5a and 8a armchair nanosheets, the critical buckling strain has a decreasing trend and the values of the critical buckling strain reduces with increasing the lengths and compared to similar zigzag nanotubes they have less buckling strain and more variation as the diameter changes.The ratio of critical buckling strain to nanotube diameter changes is always greater than one and in the diameters less than 7 angstroms, it has a decreasing trend.For the finite length armchair nanotubes obtained by wrapping 2z zigzag nanosheets, the results show that as the diameter increases the critical buckling strain decreases and the buckling strain values are lower compared to the infinite length nanotubes have greater variations with increasing diameter.Compared to armchair nanotubes obtained by wrapping similar armchair nanosheets, they had lower critical buckling strain values and lower variations.Compared to zigzag nanotubes obtained by wrapping zigzag nanosheets, they have a lower critical buckling strain but have larger diameter-dependent variations.The ratio of critical buckling strain to nanotube diameter changes is always greater than one and in the diameters smaller than 7 angstroms, it has a decreasing trend and higher than 7 angstroms, it has an increasing trend.For armchair nanotubes obtained by wrapping 3z zigzag nanosheets the critical buckling strain has a decreasing trend and its values are lower than armchair nanotubes obtained by wrapping armchair nanosheets and smaller diameter-dependent variation.Compared to zigzag nanotubes obtained by wrapping wrapped zigzag nanosheets, they have lower critical buckling strain and larger diameter-dependent changes.The ratio of critical buckling strain to nanotube diameter changes has always been greater than one and in diameters smaller than 5 angstroms, it has a decreasing trend, which is reversed in diameters higher than 5 angstroms.For armchair nanotubes obtained by wrapping 5z and 9z zigzag nanosheets, the critical buckling strain has a decreasing trend and compared to the infinite length nanotubes www.nature.com/scientificreports/they have lower critical buckling strain values, which increases with increasing length.Compared to zigzag nanotubes obtained by wrapping zigzag nanosheets, they also have lower critical buckling strain and higher diameter-dependent changes.
In Table 6, the critical buckling strain results of (n, n/2) nanotubes, which is a chirality between zigzag nanotubes with chirality of (n, 0) and armchair nanotubes with chirality of (n, n) structures, are compared with the finite and infinite lengths nanotubes obtained by wrapping armchair and zigzag nanosheets.Same as other structures, in general, with increasing diameter, the critical buckling strain has a decreasing trend.For infinite length nanotubes, the critical buckling strain of (n, n/2) nanotubes have higher trend of changes in smaller diameters.In addition, initially the critical buckling strain decreases with higher rate as the diameter increases and at larger diameters the trend of change decreases.At very small diameters the critical buckling strain diameters of (n, n/2) nanotubes is higher than zigzag and armchair nanotubes, and with increasing diameter www.nature.com/scientificreports/For (n, n/2) nanotubes obtained by wrapping 3a armchair nanosheets, the critical buckling strain has a decreasing trend with increasing the diameter and has a larger buckling strain and larger diameter-dependent changes compared to similar nanotubes with infinite lengths and compared to similar nanotubes with shorter lengths, they have lower buckling strain and diameter-dependent changes.Compared to zigzag and armchair nanotubes with similar diameters obtained by wrapping similar nanosheets, they have a lower critical buckling strain and higher diameter-dependent changes.The ratio of critical buckling strain to nanotube diameter changes is always greater than one and in the diameters lower than 7 angstroms there is a decreasing trend and in the diameters higher than 7 angstroms, it is reversed.For (n, n/2) nanotubes obtained by wrapping 5a and 8a armchair nanosheets, the critical buckling strain and diameter-dependent changes have a decreasing trend with increasing diameter and length of the structure, and compared to the infinite length nanotubes, the critical buckling strain and diameter-dependent variations are larger.Compared to similar zigzag and armchair nanotubes, they have lower critical buckling strain and larger diameter-dependent and critical buckling strain variations.By studying (n, n/2) nanotubes obtained by wrapping 2z zigzag nanosheets, it was found that the critical buckling strain decreasing trend is valid, and compared to infinite length nanotubes, they almost has a similar critical buckling strain and higher diameter-dependent variations.Compared to (n, n/2) nanotubes obtained by wrapping 2a armchair nanosheets they have lower critical buckling strain and smaller diameter dependent changes.Compared to zigzag and armchair nanotubes obtained by wrapping similar nanosheets, the critical buckling strain is lower and diameter dependent changes are higher.The ratio of critical buckling strain to nanotube diameter changes is always greater than one and in the diameters lower than 7 angstroms there is a decreasing trend and in the diameters higher than 7 angstroms, it is reversed.The critical buckling strain for (n, n/2) nanotubes obtained by wrapping 3z zigzag nanosheets has a decreasing trend and lower critical buckling strain and larger diameterdependent variations compared to infinite length nanotubes.Compared to nanotubes with smaller length, they also have lower critical buckling strain and diameter-dependent changes and Compared to nanotubes obtained by wrapping armchair nanosheets with similar length and diameter, they have lower buckling strain and diameterdependent variations.Compared to zigzag and armchair nanotubes with similar lengths and diameters wrapped in identical plates, they have smaller critical buckling strain and larger diameter dependent changes.Compared to zigzag and armchair nanotubes with similar lengths and diameters obtained by wrapping identical nanosheets, they have smaller critical buckling strain and higher diameter dependent changes.The ratio of critical buckling strain variations to nanotube diameters is always greater than one and in diameters lower than of 23 angstroms there is an increasing trend and when the diameter goes above 23 angstroms it experiences a decreasing trend.For (n, n/2) nanotubes obtained by wrapping 5z and 9z armchair nanosheets, the critical buckling strain and diameter-dependent changes have a decreasing trend with increasing diameter and an increasing trend with increasing structure length and compared to the infinite nanotubes they have lower critical buckling strain and larger diameter dependent changes.Moreover, compared to zigzag and armchair nanotubes with similar lengths and diameters obtained by wrapping identical nanosheets, they have lower critical buckling strain and larger diameter dependent changes.
The results in Tables 4, 5 and 6, also show that by increasing the length to diameter ratio, for a given chirality that leads to a particular diameter, the critical buckling strain has a decreasing trend.To compare the strength of the structure based on the length to diameter ratio, the diameter of the nanotubes must be equal because the ratio of length to diameters with dissimilar diameters has a different structural physics, that in addition to quantum effects of fine scaling, structure physics also affects the critical buckling of structure and no accurate comparison could be made.
Table 7 presents the critical buckling strain results for nanotubes with finite length and diameter, obtained by wrapping different nanosheets.In this table, in order to compare the quantum effects of finite scaling on the nanotubes with specific length, diameters are chosen in a way that covers the range of 5 to 30 angstroms, and for each specific diameter, the zigzag, armchair and nanotubes with chirality of (n, n/2) are presented.For nanotubes with a specified length obtained by wrapping armchair nanosheets, the critical buckling strain variations of the two armchair and zigzag structures are closer to each other at smaller diameters, which increase at larger diameters, and at much larger diameters where the structure of nanotubes approaches the nanosheet structure, these critical buckling strain changes decrease.The critical buckling strain variations of armchair and (n, n/2) nanotubes in smaller diameters, increase as the diameter increases, and have a decreasing trend at larger diameters.As the length of the structure increases, differences between different chiralites become more specific and there is an increasing trend.At larger lengths and diameters, the intensity of change of the critical buckling strain decreases, and this trend confirms that as the diameter increases, the strain would eventually reaches the value of nanosheets, and at larger lengths, it approaches the critical buckling strain of infinite length nanotubes.For nanotubes with specific lengths and diameters obtained by wrapping zigzag nanosheets, the results show that, same as the nanotubes obtained by wrapping armchair nanosheets, the critical buckling strain difference of nanotubes with different chirality increases with increasing length and diameter and at larger lengths and diameters, the intensity of this increasing trend decreases.By comparing the critical buckling strain of nanotubes obtained by wrapping armchair and zigzag nanosheets it was found that the critical buckling strain difference of armchair and zigzag nanotubes, in nanotubes obtained by wrapping zigzag nanosheets was higher than the ones obtained from armchair nanosheets.In addition the critical buckling strain difference of the armchair and (n, n/2) nanotubes has an increasing trend at smaller diameters and decreases at larger diameters.

Validation with other works
To show the validity of the current study, in addition to Table 3 and Fig. 13, we also represent some of the available data from other literatures in Tables 8, 9   www.nature.com/scientificreports/It should be noticed that since the effect of fine scaling has never been studied before, the available data is used to validate the none-fine scale results, which shows a good agreement.In different articles, the results show that the critical buckling force decreases with increasing the ratio of length to diameter, and the critical buckling force, which is directly related to the critical buckling strain, decreases with increasing the ratio of length to diameter.In addition, with increasing the ratio of length to diameter, the buckling strength of the structure decreases.The results and the process of its changes in this article are in good agreement with other literatures.It should be noted that for a nanotube with different lengths and diameters, the length to diameter ratio could be the same, which allows the buckling strain to be different, and for a more accurate comparison, the length to diameter ratio for structures with the same diameter should be compared to each other.Tables 8 and  9 show the comparison between exact and approximate buckling loads P cr (nN) for the simply supported based on nonlocal Euler-Bernoulli beam model and the comparison of axial buckling load of the CNT embedded in Winkler, Pasternak and Kerr's medium, respectively.In addition, e 0 a shows the effects of nonlocal parameter, L is the length and d is rod diameter.Tables 10 and 11 show the MD buckling results for SWCNTs under axial load and MD buckling results for SWCNTs with L/D ≤ 10 under torsion, respectively.In these two tables, ε cr is www.nature.com/scientificreports/ the critical buckling strain,P cr is the critical buckling load and T cr and θ cr are the Critical buckling torque and critical end rotation, respectively.

Conclusion
In this paper, the buckling behavior of carbon nanotubes with different chirality and lengths under axial loading were investigated for the first time using the quantum mechanics and molecular mechanics methods.By comparing the buckling behavior of nanotubes at different lengths and diameters, the effect of quantum effects on the buckling behavior of one-dimensional (infinite length nanotubes) and zero-dimensional (finite length nanotubes) nanostructures can be investigated.First, by combining quantum mechanics and molecular mechanics, molecular mechanics coefficients were obtained for finite and infinite length nanostructures, and then, by using molecular mechanics, the buckling behavior of finite and infinite nanotubes was investigated with respect to their diameters.In addition, besides studying the quantum effects of fine scaling dependence of nanostructures on the longitudinal changes, the quantum effects of atomic arrangement have also been investigated on the buckling behavior of the nanostructures in such a way that fine-scale nanostructures with identical lengths obtained from nanosheets with different atomic arrangement have also been studied.The results show that the critical buckling strain of the CNTs, which reflects the buckling behavior of the nanotubes, is influenced by the atomic arrangement and the type of structure, which are wrapped to make nanotubes, as well as the length of the structure.The results show that, in general, nanotubes with zigzag atomic arrangement are more resistant to axial load, which leads to buckling strain, compared to nanotubes with different chiralities and with changes in diameter and length, they show less variation.With increasing length for nanotubes obtained by wrapping armchair nanosheets, the buckling strain changes are more affected by the quantum effects than the nanostructures obtained by wrapping zigzag nanosheets.In other words, it is shown that the buckling behavior of structure with zigzag atomic arrangement is less affected by the quantum effect than the structure with armchair atomic arrangement.For nanotubes obtained by wrapping armchair nanosheets, the critical buckling strain decreases with increasing the length of the nanotubes.In other words, the results show that the smaller nanotubes can withstand higher loadings before entering the buckling process.As a conclusion, it should be noted that the buckling behavior of nanostructures at very small dimensions is highly dependent on the length, and atomic arrangement of the structure.In addition, at a very small scale, the quantum effects have an important impact on the behavior and properties of the nanostructures that if not considered in the calculations done by scientists, it would greatly cause errors in their conclusions and failure to achieve their intended purpose.

Figure 2 .
Figure 2. One-dimensional carbon nanostructures with specified width and (a) armchair and (b) zigzag boundary atoms.

Table 1 .
Effective width for investigation of mechanical properties characterizing of the fine-scale structures.

Table 2 .
Molecular mechanics coefficients for graphene structure and fine-scale nanostructures.

of structure Young's modulus (Gpa × nm) Poisson's ratio k ρ (nN/nm) C θ ( nN × nm) Cω (nNnm)
Figure 6.Buckling strain values (ε 0 ) of infinite length (n, n) nanotubes with armchair atomic arrangement with respect to inverse buckling mode parameter ξ.Figure 7. Changes of buckling strain values of armchair nanotubes with specified lengths obtained by wrapping armchair nanosheets with respect to inverse buckling mode parameter.

Table 3 .
Critical strains of single-walled CNTs as obtained by various researchers using atomistic simulations under different simulation conditions.

Table 4 .
Critical buckling strain results of zigzag nanotubes with finite and infinite lengths obtained by wrapping armchair and zigzag nanosheets.

Table 5 .
Critical buckling strain results of armchair nanotubes with finite and infinite lengths obtained by wrapping armchair and zigzag nanosheets.

Table 7 .
, 10 and 11.Critical buckling strain of nanotubes with finite length and diameter wrapped obtained by wrapping different nanosheets.

Table 8 .
Comparison between exact and approximate buckling loads P cr (nN) for the simply supported based on nonlocal Euler-Bernoulli beam model.e 0 a shows the effects of nonlocal parameter, L is the length and d is rod diameter.

Table 9 .
Comparison of axial buckling load of the CNT embedded in Winkler, Pasternak and Kerr's medium.e 0 a shows the effects of nonlocal parameter, L is the length and d is rod diameter.